Superbroadcasting of harmonic oscillators mixed states 



Giacomo M. D'Ariano,^'0 Paolo Perinotti,^'0 and Massimiliano F. Sacclii^'^'0 

' Dtpartimento di Fisica "A. Volta" and CNISM, via Bassi 6, 1-27100 Pavia, Italy. 
^CNR - Istituto Nazionale per la Fisica della Materia, Unitd di Pavia, Italy. 
(Dated: February 1, 2008) 

We consider the problem of broadcasting quantum information encoded in the displacement pa- 
rameter for an harmonic oscillator, from N to M > N copies of a thermal state. We show the 
Weyl-Ifeisenberg covariant broadcasting map that optimally reduces the thermal photon number, 
and we prove that it minimizes the noise in conjugate quadratures at the output for general in- 
put states. We find that from two input copies broadcasting is feasible, with the possibility of 
simultaneous purification {superbroadcasting). 
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On! i I- INTRODUCTION 

D , The no-cloning theorem m states that one cannot produce a number of independent physical systems prepared in 

■ identical states out of a smaller amount of systems prepared in the same state. This fact challenged the scientific 
\^ ^ community to find transformations which can approximate the cloning transformation with the highest possible 

fidelity. A whole branch in the literature is devoted to this problem, and optimal cioners have been found, for qubits 
1^ Q , for general finite-dimensional systems @ , for restricted sets of input states |^ Q , and for infinite-dimensional 
^ . systems such as harmonic oscillators — the so called continuous variables cioners . In all these cases pure states have 
been considered. On the other hand, the problem of cloning mixed states is somehow more interesting, since a more 
general type of cloning transformation can be considered — the so-called broadcasting — in which the output copies are 
in a globally correlated state whose local "reduced" states are identical to the input states. This possibility has been 
[ considered in Ref. , where it has been shown that broadcasting a single copy from a noncommuting set of density 
\^ ■ matrices is always impossible. The no broadcasting theorem was the only interesting result for mixed states for many 
years, and this led to the belief that a generalization of the no-cloning theorem to mixed states should hold. However, 
more recently, for qubits an effect called superbroadcasting jlQI] has been discovered, which consists in the possibility 
of broadcasting the state while even increasing the purity of the local state, for at least > 4 input copies, and for 
^ , sufficiently short input Bloch vector (and even for A^ = 3 input copies for phase-covariant broadcasting instead of 
' universal covariance and these superbroadcasting maps can be mixed with depolarizing maps in order to get 

Ci [ exact broadcasting. 

^ ■ In the present paper, we analyze the broadcasting of mixed states of an harmonic oscillator by a signal-preserving 
map. More precisely, this means that we consider a set of input states obtained by displacing a fixed mixed state by 
a complex amplitude in the harmonic oscillator phase space, while the broadcasting map is covariant with respect to 
the (Weyl-Heisenberg) group of complex displacements. We will focus first on displaced thermal states (which are 
^ , equivalent to coherent states that have suffered Gaussian noise), and then we will show that all results hold for any 

■ covariant set of mixed states in terms of noise of conjugated quadratures. 

As we will see, superbroadcasting is possible for harmonic oscillator mixed states |l2l Il3l |. namely one can produce 
a larger number of copies, which are locally purified on each output oscillator, and with the same signal as the input. 
For displaced thermal states, for example, superbroadcasting can be achieved for at least N = 2 input copies, with 
thermal photon number riin > ^, whereas, for sufficiently large ni„ at the input, one can broadcast to an unbounded 
number M of output copies. For purification (i.e. M < A), quite surprisingly the purification rate is rioutlnin = A~^, 
independently on M. The particular case of 2 to 1 for noisy coherent states has been reported in Ref. We will 

prove also that perfect broadcasting is possible, provided that one knows the input thermal photon number. 
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II. COVARIANT BROADCASTING FOR THE WEYL-HEISENBERG GROUP 

We consider the problem of broadcasting N input copies of displaced (generally) mixed states of harmonic os- 
cillators (with annihilation operators denoted by ag, ai, a^v-i) to M output copies (with annihilation operators 
bo, bi, ^M-i)- In order to preserve the signal, the broadcasting map ^ must be covariant, i. e. in formula 

.^{D{af^ED{ay^^) = D{a)'^^' 3§{E)D(a)^®^^ , (1) 

where Dc{a) = exp(ac^ — a*c) denotes the displacement operator, and S represents an arbitrary A^-partite state. It 
is useful to consider the Choi-Jamiolkowski bijective correspondence of completely positive (CP) maps 3§ from Tijn 
to Tiout and positive operators Rgg acting on Hout ® "Win , which is given by the following expressions 

R^ = ,'^(^^{\n){n\) , 

^(p) = Tri„[(/out®p")i?.^] , ^ ' 

where jfi) = X^^^o \i^ri)\''Pn) is a maximally entangled vector of Ti-®^ , and X'^ denotes transposition of X in the basis 
IV'n)- In terms of the operator Rag the covariance property ^ can be written as 

[Rss,D{a)®^^' ®D{a*)®^]^Q , Va e C . (3) 

In order to deal with this constraint we introduce the multisplitter operators Ua and J7f,, that perform the unitary 
transformations 

N-l 



^ 1=0 
M-l 

M 



M-l 

U,b,ul - ^ E ^"^b^ ■ (4) 



(=0 



Notice that such transformations perform a Fourier transform over all input and output oscillators. Moreover, we will 
make use of the squeezing transformation Saoba defined as follows 

[Saobo^Cln] = [Saobo^K] =0, n > 

Saob„alsl^f,^ ^ fial - i^bo , (5) 
SaobobaSi^^^^ = /i6o - i^al , 



with /i = AI/{M — N) and ly = ^ N/ (M — N). The squeezing transformation here acts as an hyperbolic transfor- 
mation for just oscillators aq and feoi by leaving all other oscillators unaffected. In terms of such operators, condition 
(0) becomes 

[SlbM ® UDRp^Ub ® Ua)Saob„,Db„{VM^a)] = . (6) 

Hence, upon introducing an operator B on oscillators bi, 5a/-i, clq, cln-Ii the operator R^gg can be written in the 
form 

R.m = (C/fc ® U,)Saobo {ho ® B)Sl,^ {Ul ® C/t). (7) 

Notice that i?,^ > is equivalent to i? > 0. The further condition that ^ is trace-preserving in terms of R^ 
becomes Trh[i?^] = la, b and a collectively denoting all output and input oscillators, respectively. From the trace 
and completeness relations for the set of displacement operators, namely / (Pa D{a)AD^ (a) ~ Tr[^]/, and A — 
J (PaTT[D^{a)A]D{a), (see, e.g., Ref. [13), the condition Trb[i?,®] = la is verified iff 

M-l \ /M-l \ 

Db. m Saobo {ho ® B)Sl. {(^Dlm]^I. (8) 




From the relation Dbo{Po)Saobo = SaoboDbo{fJ-Po) ® D\^{v(3q), one obtains the condition 

^^b/bo,ao[B]^v'^Ia/ao , (9) 
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where ajai denote all the input oscillators apart from a^, and similarly for hjhi. 
We will now consider the map corresponding to 

i? = i^'|0)(0|6/,„® |o)(ou„®4/,„ . (10) 

Applying the corresponding map to a generic A^-partite state S we get 

^(S) - Tr,[(/, ® E^){Ut ® Ua)Sa„bo{ho ® B)Sl,^^{Ul ® [/])] , (11) 

which is equivalent to 

^(S) = TrJ(/, ® UlE^Ua){U, ® Ia)Saobo {ho ® S)^I„6„ (C/fe^ ® /a)] • (12) 

Using the expression in Eq. (|10(l we obtain 

^(S) = [/fc {Tr,J(/,„ ® CJ^aofcoa^o ® |0)(OUo)^l„,J ® |0)(0|fc/b„} C/,^ , (13) 
where ^ = Tra/ajC/^S^C/a]. Notice that 



e = / ^i?(7)^Tr[(i?,„(7)t®4/.jC/tS-L/J 
^i?(7)^Tr[C/:(I?„,(7)* ® /a/ao)C/:S], 



Pa 



n + 1 \n + 1 



(14) 



(15) 



and using the complex conjugate of Eq. Q we have 

—^D{jYTt[D{j*/VN)^''E] 

TT 

Now, we can easily evaluate Saoboiho ^ |0)(0|ao)'S'aof,oJ by expanding the vacuum state as 

|0)(0|ao= / ^e-'^DM, (16) 

obtaining 

Saoboiho ® |0)(0UJ5l^_,^ = 1 ^^e-^Ao(^7*) ® ^ao(/^7) ■ (17) 

Hence, Eq. (|13|l can be rewritten as 

^f/,(A„(7*) ® |0)(0|,/,jC/te-^Tr[i?,„(Ai7/^)r] • (18) 
As an example, we will now consider A displaced thermal states 

• ^ Dia)(^X'^ D{a)\ (19) 



from which we want to obtain M states, the purest as possible. Thanks to the covariance property, it is sufficient to 
focus attention on the output of Pq^ ■ For a tensor product of thermal input states S = Pq^ , exploiting the fact that 

Uli'E'j'jo^ a]aj)Ua = 'EJJo^ a^Qj, we have 

e = r = po , (20) 

and recalling the following expression for the thermal states 



1 / n 



n + 1 \ n + 1 



""'le~^('"+i)L'(/3), (21) 
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we obtain 



-[m (2n+l) + l] 



7Ubih){l\bo®\0){0\b/bo)Uy 



(22) 



f|7/VM)(7/VM|®^^e 



/l7)(7l 



where 



1.9. X 1 2Mn + 2M - N 

2n' + 1 = — [^2(2n + 1) + l] = — 



..... .r ■ (23) 

The above state is permutation-invariant and separable, with thermal local state at each oscillator with average 
thermal photon 

_„ _ n' _ Mn + M - N 
^ ^M^ MiV ■ 
More generally, for any state S, the choice pU|) gives M identical clones whose state can be written as 



(24) 



p'^ J lii!e-^(w-*+i){Tr[Si?t(a/Ar)®A']}i^(a) . (25) 
Since for any oscillator c one has 

Ax^ + Ay,^ = i + (cM-|(c)p, (26) 

it is easy to verify that the superbroadcasting condition (output total noise in conjugate quadratures smaller than the 
input one), is equivalent to require a smaller number of thermal photons at the output than at the input, namely 



Mn+M-N _ M-N 

n > ^ n> 



MN 

This can be true for any > 1, and to any Af < 00, since 

M - N 

lim 



M{N - 1) 



M^oo M{N - 1) N-1 



> 0. 



(27) 



(28) 



III. PROOF OF OPTIMALITY FOR THE CHANNEL IN EQ. JTTll 



Actually, the solution given in Eq. H24|l is optimal. To prove this, in the following we will show that the expectation 
of the total number of photons Tr[^;^g ^ b]bi^{pf^)] of the M clones of p cannot be smaller than Mn". Since the 
multisplitter preserves the total number of photons we have to consider the trace 

= Tr h\bi ® {UIpI^'^Uo) Sa^boiho ® B)Sl,^ . (29) 



1=0 



We can write W = Wq + J2fii ^ with 



Wq = Tr 
Wi = Tr 



Sl,, ({hlho ® Ib/bo) ® (UU^'^Ua)) SaMho ® B) 

Slbo (ih/b, ® blbi) ® ([/tp5«^C/,)) SaMho ® B) 



(30) 



for 1 < / < M — 1. Now, since Wi > Q, W > Wq. Moreover, using the identity c^c ~ ~daa*e 2 Dc{a)\a=a'=o, one 
obtains 



Trfco Sl,bo (^0^0 ® (j) Saobo 



— ~B . 



TT 



Trfc„ [Dbo (pa Da„ (pi - vo*)\ Tr[i:)(7) V]e^ U=a*=o 

_ ajao + ^^Tr[aJaofT] + 1 



(31) 



OL — Oi* 
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then, from Eq. (j^l and positivity of B, one has 

_Tr[(/,/b„ (g> alao (g> Ia/a„){h/bo ® (Uapf^UlrjB] 



Wo 



^i'^T[{Ib/bo ^ lap ^ TTa„[alao{Uapi^Ulr])B] + ^2 (32) 



n^n + l N M - N 
> - — —n H Mn". 

- M N 



In fact, one can easily check that the choice of B in Eq. (|10|) saturates the bound (|32|l . Notice that for n = one has 

\I-N 
MN ' 



N coherent states at the input, and n" = ^jri^, namely one recovers the optimal cloning for coherent states of Ref. 



Also the more general solution given in Eq. I|25() for arbitrary input state is optimal, in the sense that it represents 
the state of M identical clones with minimal photon number, which is given by 

rr \hU n Tr[atapo] 1 1 

Tr[6t6H = ^^ + ^-^. (33) 



IV. EXACT BROADCASTING 

In the previous section we proved that it is possible to superbroadcast thermal states provided that fi > 
Here we will prove that in this case it is also possible to broadcast perfectly, namely with local fidelity equal to 1. 
Suppose indeed that the Choi-Jamiolkowski operator is given by Rag = (Ut (8) Ua)Saobo {h ® B')S\^i^^{Ul ® Ul), with 

B' = V'' [pT^^),I,, ® |0)(0|ao ® la/ao- (34) 

Then by using Eq. (|21|l we replace the first line of Eq. (|22|) by 



IT" 

a^(2ft' + l)- ^'° = i3 (2m+l) 



X e 



(35) 



X e 

,A/-1 



^ ^'1 (2fl'-2m)- ^fc=",'^''' (2m+l) 



where 7fe = J2j=o ^ " ■> ^^'^ Jacobian for the change of variables is clearly 1 . Taking the partial trace 
over all output spaces but one, we get the following local state on the fc-th oscillator 

^'^ /^^(^)e-^(2-+i+^), (36) 



which is a thermal state with thermal photon number 

„ n' + miM-l) n 1 1 M-1 
n" = — - — i '- = \ h m . (37) 

M N N M M ^ ' 

In order to have n" = n it is sufficient to choose a suitable thermal number to, which is the solution of the following 
equation 

n 1 1 Af - 1 

n = \ hTO , (38) 

N N M M ' ^ ^ 



namely 



_ _ M{N-l)n _ (M-N) 

7V(Af-l) iV(Af-l)' ^ ' 
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which is positive iff 



M-N 

This definitely proves that perfect broadcasting is possible provided that superbroadcasting is. On the other hand, 
since superbroadcasting provides the maximum output purity, it is impossible to have perfect broadcasting when 
n < jvjp^zY): namely either superbroadcasting and perfect broadcasting are both possible, or they are both impossible. 

V. NOISE IN CONJUGATE QUADRATURES 

We will consider now an alternative derivation, which involves radiation modes and makes use of a theorem for 
linear amphfiers. We are interested in a transformation that provides M (generally correlated) modes 60, bi, 6m-i 
from N uncorrelated modes ao,ai, ...,aAr-i, such that the unknown complex amplitude is preserved and the output 
has minimal phase-insensitive noise. In formula, we have input uncorrelated modes 

(flj) = a , 

Axl + Ayl ^ J, > ^, (41) 

for all i = 0, 1, .., — 1, where Heisenberg uncertainty relation is taken into account. The output modes should satisfy 

(bi) = a , 

Axl + Axl - r > i , (42) 

and we look for the minimal T. The minimal T can be obtained by applying a fundamental theorem for phase- 
insensitive linear amplifiers 17] : the sum of the uncertainties of conjugated quadratures of a phase- insensitive amplified 
oscillator with (power) gain G is bounded as follows. 

AXl + AYi > G{AXl + AYl) + ^ , (43) 

where A and B denotes the input and the amplified mode, respectively. Our transformation can be seen as a phase- 
insensitive amplification from the mode A = -^^^^^ cli to the mode B = "^;^X]f=o^^» with gain G = and 
hence Eq. H43|) should hold. Notice that generally for any mode c one has 

Ax^-f Ay2 = i + (ctc>-|(c)p. (44) 

Hence, the bound can be rewritten as 

{B^B) - > G{{A^A) + 1 - - 1. (45) 

In the present case, since modes are uncorrelated, one has 

JV-l (n-1 

i j'=0 y 1=0 i^j 

= (7 + lap -\) + {N- 1)|«P = 7 + - \. (46) 

where 7 = -i- ^^q^ 7;, and so the bound Eq. 143|) is written as 

(StS) >G(7+i)-l + M|a|2. (47) 

On the other hand, one has 

1 1 I 1 

(^^^) - M E (^^^^■) ^ M E ^^b\br)^bh) - ^^(r + - (48) 
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Eqs. H47|l and (|48|l together give the bound for the minimal noise T 

T--> —h- -) + — -— . (49) 

The example in the previous sections corresponds to 7 = n + i and T — n" + ^. A similar derivation gives a bound 
for purification, where N > M. In such a case G < 1, and Eq. 143() is replaced with 

AXl + AFj > GiAXl + Ayj) + , (50) 

and one obtains the bound 

r-i>>-i). (51) 

We would like to stress that the derivation of all bounds in the present section relies on the theorem of the added noise 
in linear amplifiers, namely only linear transformations of modes are considered. Hence, in principle, these bounds 
might be violated by more exotic and nonlinear transformations. Therefore, the derivation of Eq. (|32|) is stronger, 
since it has general validity. 

By a similar derivation, using the bound for phase-conjugated amplifiers AXj^ + AYj > G{AXj^ + AY^) + ^^y^, 
one can obtain the bound for phase- conjugation broadcasting 

r4>^(7 + i)^ (52) 



VI. EXPERIMENTAL IMPLEMENTATION 



The optimal broadcasting can be easily implemented on radiation modes, by means of an inverse A^-splitter which 
concentrates the signal in one mode and discards the other — 1 modes. The mode is then amplified by a phase- 
insensitive amplifier with power gain G = Finally, the amplified mode is distributed by mixing it in an M-splitter 
with M—1 vacuum modes. Each mode is then found in the state of Eq. (|25|l . In the concentration stage the A^ modes 
with amplitude (a;) = a and noise Axf + Ayf = ji are reduced to a single mode with amplitude V^TVa and noise 7. 



The amplification stage gives a mode with amplitude \fMa and noise 7' = 7^ + ^ — \- Finally, the distribution 
stage gives M modes, with amplitude a and noise T = (7' + ^'^2~^ ) ^ach. In Fig. ^we sketch the scheme for 2 to 
3 super broadcasting. 

In Ref. |l8l | it was shown experimentally that phase insensitive amplification can be obtained by a setup consisting 
of a beam-splitter, a heterodyne detector and a conditional displacement. In the following we give an algebraic 

derivation of this result. Consider a mode in a state p — J ~~/(7)£'(7) coupled to another mode in the vacuum 
through a beam-splitter with transmissivity r. The output is given by the bipartite state a 



(53) 



where we performed the change of variables p ^ t/3 + Vl — t^7, 7 ^ r7 — Vl — Now, the reflected mode 

is measured by heterodyne detection, and conditionally on the measurement outcome a, a displacement D(ka) is 
performed on the transmitted mode, whose state is then given by 



P' = 



TT 

d^ad-pd-j ^_ ^.,-V^., /(r7+ VT^/3)e"(^-*-''*)--^-i?(7)e- 

<PP<Pl .{2)1, a\ 



e- ^ ' " /(r7 + - T^(3)D{ka)D{j)D{ka)^0\D{ayD{f3)D{a)\0) 



S'-^\kj-P)e- -^'^ /(t7+ Vl-T^/3)e'^ = 

71 - 

^ /(7(r + A;v/T^))e-4^[^'+(^— ^^)'li?(7). (54) 
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On the other hand, the action of a phase-insensitive amphfier on p can be easily calculated and produces the partial 
output state 

/ = / ^e-^/(A*7)i^(7). (55) 

The following conditions 

^ = r + fcVl-r2, v"^ = k^ + (fcr- y/l-T^f, fi^ - ly^ = 1, (56) 

which equivalent to 

k ~ V, T = —, (57) 

imply that p' — p". Hence, by tuning the beam splitter transmissivity and the parameter of the conditional displacement 
fc, one can then simulate the amplifier by a linear device assisted by heterodyne and feed-forward. 




FIG. 1: Experimental scheme to achieve optimal superbroadcasting from 2 to 3 copies. This setup involves just a beam splitter, 
a phase-insensitive amplifier and a tritter, which in turn can be implemented by two suitably balanced beam splitters. The 
phase-insensitive amplifier can be implemented by a beam splitter and heterodyne-assisted feed-forward. The output copies 
carrying the same signal as the input ones are locally more pure, the noise being shifted to classical correlations between them. 



VII. CONCLUSION 

In conclusion, we proved that broadcasting of M copies of a mixed radiation state starting from N < M copies 
is possible, even with lowering the total noise in conjugate quadratures. Since the noise cannot be removed without 
violating the quantum data processing theorem, the price to pay for having higher purity at the output is that the 
output copies are correlated. Essentially noise is moved from local states to their correlations, and our superbroad- 
casting channel does this optimally. We obtained similar results also for purification (i.e. M < N), along with the 
case of simultaneous broadcasting and phase-conjugation, with the output copies carrying a signal which is complex- 
conjugated of the input one. Despite the role that correlations play in this effect, no entanglement is present in the 
output (as long as the single input copy has a positive P- function), as it can be seen by the analytical expression 
of the output states. Moreover, a practical and very simple scheme for experimental achievement of the maps has 
been shown, involving mainly passive media and only one parametric amplifier. The superbroadcasting effect has 
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a relevance form the fundamental point of view, opening new perspectives in the understanding of correlations and 
their interplay with noise, but may be also promising from a practical point of view, for communication tasks in the 
presence of noise. 
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